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We discuss an extended standard model electroweak sector which contains a stable 



scalar dark matter particle, the D boson. To search for the D boson at the LHC 
we exploit the flavor-changing neutral current (FCNC) top quark decay, t — > cDD, 
mediated by the lightest standard model-like Higgs h in a two Higgs doublet model 
framework. The branching ratio for t — > cDD in this case can be as high as 10 3 , 
^■i, after taking into account constraints arising from the D boson relic abundance. With 

an integrated luminosity of 10 (100) fb _1 , the 14 TeV LHC can explore values of 
this branching ratio that are one (two) order of magnitude smaller in tt production 
^ ! with tt -> cU~(cbl + ) +Xt- For a D boson mass < 60 GeV, m h o < 2M Z , 10 fbr 1 

| luminosity and a branching ratio BR(t — > cDD) ~ 10~ 4 , the estimated number of 

signal events at the 14 TeV LHC is of order 80. 



I. INTRODUCTION 

•rH 

X 

A large number of direct and indirect experiments are currently underway |2Nl3| search- 



ing for the weakly interacting massive particle (WIMP) whose relic abundance presumably 
provides about 23% of the universe's energy density l|. The successful launch of the Large 
Hadron Collider (LHC) at CERN provides an unparalleled opportunity to produce WIMPs 
in p-p collisions and infer their existence through large missing energy events. The interplay 
between experiments at the LHC and the direct and indirect searches will play a crucial role 
in identifying the WIMP dark matter particle. 

It is almost universally agreed that the Standard Model (SM) offers no viable WIMP 
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candidate, and therefore some extension of this highly successful theory is warranted. One 
particularly simple extension is to add a SM singlet real scalar field which yields a spin zero 
particle with mass on the order of the electroweak scale or less Jjj]. An unbroken Z2 
parity, under which only the scalar field is odd, makes this spin zero particle (called D boson 
here) stable. For recent discussions see 16M24 1 . 

At the renormalizable level D only couples to the SM Higgs doublet. This coupling must 
be carefully adjusted to reproduce the required relic density of D, while making sure that 
constraints arising from the direct searches are not violated. However it is hard to achieve 



this within the SM+D framework 



21 



221 ] . In order to obtain a consistent scenario with D 



boson dark matter, it is desirable to consider an extension of the SM, such as the two Higgs 
doublet model (2HDM) that we discuss here 21]. 

In this paper we propose a search for the D boson at the LHC by considering the impact 
D could have on rare top decays. With a total cross section a(tt) ~ 800 pb at LHC, a 
large number of ti pairs will be produced and top quark physics will be studied in great 
detail. In particular, flavor changing neutral current (FCNC) decays of the top quark such 
as t — > c/i°(7, Z, g), with branching fractions as low as 10~ 5 or so, can be explored j^. In 
the presence of D, one could envisage FCNC processes such as t — > ch° — > cDD which, 
unless highly suppressed, should be taken into consideration. In the SM+D model, such 
processes arise at the loop level and are heavily suppressed. We therefore consider as a 
concrete example a 2HDM+D model in which the FCNC process t — > cDD arises at tree 
level, mediated by the lightest SM-like Higgs boson h°. AD boson with mass ^ 100 GeV 



in this mode 
searches 21 



is a plausible dark matter candidate which is compatible with the direct 



22] . With the parameters of the model rather tightly constrained in order 



to achieve this, the 2HDM+D model, as we will show, gives rise to some rather unique 
signatures arising from t — > cDD which may be detected at the LHC. 

The paper is organized as follows. In Section II we describe the 2HDM+D model. The 
constraints from the relic abundance of D and the FCNC top decay into a pair of D's are 
discussed in Section III. The prospects of discovering the signal associated with this process 
at the LHC are outlined in Section IV. Our findings are summarized in Section V. 
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II. SCALAR DARK MATTER IN TWO HIGGS DOUBLET MODEL 



The renormalizable interaction of a real scalar dark matter particle D boson with two 



Higgs doublet fields H\ , H 2 can be written as 23 ] 



- C D = ~^-D A + ^D 2 + D 2 (X 1 H\Hi + \ 2 HlH 2 + HH\H 2 + H^)). (1) 

Note that an unbroken Z 2 symmetry under which D — > —D has been imposed to keep the 
D boson stable. Since D couples at the renormalizable level only to the Higgs doublets, it 
interacts weakly with the rest of SM fields and plays the role of stable WIMP dark matter. 
The two Higgs doublets, after electroweak symmetry breaking, have physical components 
= (— sin (3H + , [v\ + cos aH — sin ah — ism f3A)/\/2) and = (cos (3H + , (t> 2 + sin aH + 
cos ah° + i cos f3A)/y/2). Here tan/3 = v 2 /v 1 is the ratio of the vevs of the two Higgs doublets 
and a is the mixing angle of the CP-even neutral Higgs fiels. With Z 2 unbroken, the D 
particles can only be produced or annihilated in pairs through Higgs exchange. Using the 
above information, we obtain the mass of D and the h°DD interaction (note that h° is the 
SM-like Higgs in our discussion), 

m 2 D = ml + v 2 (\i cos 2 f3 + A 2 sin 2 /3 + 2A 3 cos /3 sin (2) 
—Ch°DD = [—Xicos (3 sina + X 2 sm (3 cosa + \3Cos((3 + a)]vh°DD = \hvh°DD. (3) 

Here v 2 = v\ + v 2 = (246 GeV) 2 , and both mo and the effective coupling are free 
parameters in this model. The couplings of H and A to D are: —Chdd — (Ai cos f3 cos a + 
A2 sin /3 sin a + A3 sin(/3 + a))vHDD = XhvHDD and —Cadd = 0. For concreteness, in our 
numerical analysis we will neglect any contributions from H either by requiring a sufficiently 
small Xh or an appropriately heavy mass m#. 

A two Higgs doublets extension of the SM is denoted as 2HDM I, 2HDM II and 2HDM 
III, where 2HDM I means that only one linear combination of H\ and H 2 provides masses 
to both up and down type quarks. In 2HDM II H\ provides masses both to down type 
quarks and charged leptons, and H 2 to the up quarks. Finally, in 2HDM III, both Hi and 
H 2 provide masses to up and down type quarks, and charged leptons. In 2HDM I and II, the 
FCNC effects are generated at one loop level, and hence the FCNC top quark decay rate is 
too small to be detected at hadron colliders, even though it can be substantially larger than 
that predicted by the SM. In contrast, 2HDM III offers the possibility of a large detectable 
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rate because of the presence of tree level FCNC We therefore only focus on 2HDM III here, 
and we will refer to this model as 2HDM III+D. 



The Yukawa couplings of Hi, if 2 to the fermions in this model are given by 28] 



— C-iii = QlK^iUr + Ql^HzUr + Ql^H 1 D r + Q L X d H 2 D R 

+ UX[HiE R + T L \\H 2 E R + h.c, (4) 

where Hi = ia 2 H*. The coupling of the SM-like Higgs h° to fermions reads 

- Cu, = U L M"U R ^Lh» - VJfU^=£ih' - D L M'D R ^-h" (5) 
v sin p v sin p v cos p 

— ~ ■ ^ cos(a — B) , n — , ,, ,_, sin a , n — cosfa — /3) , n 
+ D L M d D R ^ -f±h° - E L M l E R -h° + E L M l E R ^ -f±h° + h.c. , 

v cos p v cos p v cos p 

where M u ' d ' 1 = (A"' d ' z wi + X^ ,d,l V2)/ \/2 denote the diagonalized masses of the up and down 
type quarks and charged leptons. The off-diagonal entries M u = \™v/\/2 and M d ' 1 = 
are not fixed. 

In our discussion, we follow Ref. 29| and parameterize the off-diagonal entries to have the 
geometric mean form M".' ' l = p^f' 1 ^rrii^j with p^ ~ 1 for concreteness, and pa negligibly 
small for illustration. With this parametrization for M l , the Yukawa couplings are identical 
to those in MSSM, if the off-diagonal elements are set equal to zero. To simplify our analyses, 
we further ignore the off-diagonal elements except those involving the top quark. This 
parametrization, together with the assumption of a sufficiently heavy non SM-like Higgs H 
allows one to satisfy a variety of experimental constraints, for instance from quark flavor 



changing processes and rare B decays |30|. Note that the couplings of h° to W, Z in 2HDM 
III is given by 

2M 2 M 2 

C h0ww = K sin (£ _ a )h°W 2 , C h o zz = sin(0 - a)h°Z 2 , (6) 

v v 

which will alter the Higgs decay width from its SM value. 



III. DARK MATTER CONSTRAINTS AND TOP QUARK FCNC DECAY IN 

2HDM+D MODEL 

The annihilation of a pair of -D's into SM particles proceeds through s-channel h° ex- 
change. Let us first consider DD — > h° — > //'. We parameterize the Higgs-fermion and 
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Higgs-.D interactions as 

-C Y = a{Jlf R h + h.c. + bh°D\ f = u,d,l (7) 
where R(L) = (1 ± y B )/2. In the 2HDM III+D we have 



a «. = «. ™ _ jfijSE!^, (8) 

13 lJ vsmf3 3 vsm/3 

4 = -Mi^ + Mi^^l, (9) 
3 v cos /3 3 v cos /3 



sin a ~ ; cos(a — 8) 

a, = -Ml a+ M i " IT 1 ! 10 

13 v cos 8 13 v cos 8 

X h v. (11) 



The partial decay width of h° into fermions is given by 

r(h° -+ //') = hY, N >ff\ 2 ^ - ^)f ,2 ^r + 



-3- ^ A^^|a//'| 2 (m^o — m 2 — m 2 , — 2mfmp)J (m 2 — m 2 — m 2 ,) 2 — 2m 2 m 2 ,}, (12) 



where iVj? is the number of colors of the f-fermion (3 for a quark and 1 for a lepton). The sum 
is over fermions with rrif < mo- In the non-relativistic limit the total averaged annihilation 
rate of a DD pair is then given by 

8^ r(fr° -> X') 

(<Wrel) - VannVrel ~ (4^ _ ^2 + m 2 f 2 ^ > ( 13 ) 

where r(/i° — > X') = X^E(/i° — > Xi), with /i° being a "virtual" Higgs having the same 
couplings as the Higgs h° to other states, but with a mass of 2m p. The Xi indicates any 
possible decay mode of h°. Note that the sum should also include other decay channels, 
for instance h° — > 77, gg and h° — >■ W + W~, ZZ, if allowed by the relevant kinematics. For 
a given model, T(h° — > X') is obtained by calculating the h° width and setting the mass 
equal to 2m d- In Eq. (I13jl . v re \ is the average relative velocity of the two D particles. For 
cold dark matter the velocity is small, and therefore to a good approximation, the average 
relative speed of the two D's is v re i = 2p Dcm /mo ) and s — (p/ + pf) 2 is equal to Am 2 D . 

The present relic density of D is given by = m^Sgl^, where s = 2889.2 cm" 3 is 
the present entropy density. is the asymptotic value of the ratio tid/so, with Y^ 1 = 
0.26Ay/g^Mp l m D (a ann v re i)xJ 1 through the time (temperature) evolution which is obtained 
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by solving the Boltzmann equation, where Xj = mn/Tf and Tj is the freeze-out temperature 
of the relic particle. The relic density can be expressed in terms of the critical density 

1.07 x 10 9 GeV" 1 x f 1 



Q D h 



M 



pi 



iVrel) 



(14) 



where g* is the number of relativistic degrees of freedom with mass less than Tj. The 
freeze-out temperature Xj- can be estimated through the iterative solution of the Boltzman 



equation 



33] 



In 



c(c + 2) 



45 g M Pl m D {a ann v re i) 
8 2tt 3 




In 



0.038Mpim D (a ann v re i) 



vg^J 



(15) 



\jg* x f 

where the constant c, of order unity, is determined by matching the late-time and early-time 
solutions, and g is the weak interaction gauge coupling constant. 

It is important to note that in the SM+D model, a D boson mass range 10 GeV < < 
(50, 70) GeV, with a SM Higgs mass of (120, 200) GeV, is ruled out by the upper limits on 
the WIMP-nucleon spin-independent elastic cross-section from the XENON10 and CDMSII 



experiments 
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22j . However, it has been shown that the direct detection constraints can 



be evaded if the Higgs-nucleon coupling happens to be sufficiently small due to cancelations 
among the various contributions arising from the underlying Yukawa couplings. This can 
be realized in the 2HDM+D model by setting tan a tan /3 ~ 0.405, without violating the 
relic density constraint 2l|. It is shown in Ref. 2l| that by setting tan a tan (3 = 0.45 
in the 2HDM+D model, and for all reasonable values of the D boson mass, the D boson- 
nucleon elastic cross section can be smaller than (9(10 -44 ) cm 2 , which is the upper limit from 
XENON10 Q and CDMS II 3, \ For m D > 40 GeV and with m h o > 120 GeV, the relevant 
cross section could be smaller than (9(10~ 45 ) cm 2 , which is the projected sensitivity of 
XENON100 and SuperCDMS 32]. Thus, we employ tan a tan p = 0.45 in the following 
analyses. 

For given values of and Q^h 2 , Xf and g* and therefore also (c ann v re i), can be deter- 
mined. One can then estimate the interaction strength b in Eq. ([7]). In Fig. [T]we plot Xf 
(left panel) and (cr ann v rfi ]) (right panel) versus m^, with 0.108 < fi^/i 2 < 0.1158 from cos- 
mological observations 1|. In Fig. EJ we show the allowed range for the parameter Xh = b/v 
as a function of rri£> for several values of Higgs mass m^o, and with tan /3 = 3 and 30. The D 
boson mass, we note, can be as low as 1 GeV or so. Since we are interested in producing D 
particles and studying their properties through the FCNC top quark decay at the LHC, we 
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FIG. 1: Xf (left) and (a ann v re i} (right) vs. D mass mn with 0.108 < Oo/i 2 < 



1158 fl. 



limit ourselves to a -D mass below 100 GeV. Note that as the D mass decreases, A^ becomes 
larger. For small enough mo, A^ can approach unity, which may spoil the applicability of 
perturbative calculation. Thus, we will only consider A^ < 1. 

We next explore D-physics through the FCNC decay of top quark, where a major differ- 
ence between 2HDM III+D and SM+D can show up. The decay amplitude for fa — > fjDD 
is given by 

2b 



M(ft -> fjDD) 



, —JAaiR ■ a{jL)j). 
mi + iT h m h a J J 



(16) 



In the SM the branching ratio BR[t — » cDD) was estimated to be < 10~ 13 in Ref. [23|. 
Using Eq. ([7]), the corresponding results for the 2HDM III+D model are shown in Fig. |3j 
We find that the branching ratio BR(t — > cDD) for this case can be as large as ~ 10 -3 , if 
tan /3 is sufficiently small tan /3 = 3 and the h° mass is below the h° — > VV threshold (V 
stands for vector bosons W and Z). With tan/3 = 30, the upper limit for BR(t — > cDD) 
is ~ 10~ 5 because top quark FCNC coupling decreases for larger tan/3 values. If the h° 
mass is larger than the VV threshold, we find BR{t — > cDD) < 10 -5 for small tan /3, and 
BR(t ->■ cDD) < 10" 7 for large tan/3. 



IV. OBSERVABILITY OF FCNC TOP DECAY t -»• cDD AT THE LHC 

In the following we discuss the search for D particles through FCNC top decay at the 
LHC. We are interested in the ti pair production pp — > tiX, with one of the top quarks 
decaying into a pair of D bosons through the FCNC process t — > cDD (or i — > cDD). 
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tairar taq£=0.45, tan/?=3 



tana tzq0=O.45, tan/?=30 
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FIG. 2: A fe in 2HDM III+D model vs. m D with tan a tan /3 = 0.45, tan/3 = 3 (left) and tan/3 = 30 
(right), where the shaded areas 1 (black) and 2 (red) are for m^o =120 and 200 GeV, respectively 
(same for other figures). In 2HDM III+D, we have assumed the physical Higgs h° to be much 
lighter than the other neutral scalar bosons. 
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FIG. 3: The branching ratios of £ — > cDD in 2HDM III+D as a function of m£> with tan a tan /3 = 
0.45, tan/3 = 3 (left) and tan/3 = 30 (right). T t denotes the total decay width of top quark, 
dominated by t — > bW. 

To circumvent potentially large QCD backgrounds, we require that the W boson from the 
second top quark decays leptonically. Consequently the process we are interested in is 



tt-tcb i~(cb i + ) + Xt, £ = e,[i. 



(17) 
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The overall branching fraction is given by 



BR(ti -»• r 6 c (£+ 6 c) +Xt) = 2 x - x BRt^ C DD x (1 - BR t ^ cDD ) 



(18) 



where the factor | is the leptonic decay branching ratio of the boson. 

For our numerical analyses, we adopt the CTEQ6L1 parton distribution function and 
evaluate the SM backgrounds by using the automatic package Madgraph. We work at the 
parton-level, but simulate the detector effects by the kinematical acceptance and employ 



Gaussian smearing for the electromagnetic and 
basic acceptance cuts for the event selection 



radronic energies. We employ the following 



34 



35| 



pr(£) > 15 GeV, \r](£)\ < 2.5, 
Pt(j') > 25 GeV, \ V (j)\ < 3.0, 



ARjj, AR je > 0.4, 



> 30 GeV. 



(19) 
(20) 
(21) 
(22) 



To simulate the detector effects on the energy-momentum measurements, we smear the 
electromagnetic and jet energies by a Gaussian distribution whose width is parameterized 



as 



35] 



deal 



E e /GeV 

O'had 



hah O-cal = 10%, b ca l = 0.7%, 
bhad, 0>had = 50%, bhad — 3%. 



(23) 
(24) 



V^TGeV 

In principle, the leading SM background to our signal is from the decay of W to lepton 
plus two jets. For instance, the leading irreducible backgrounds to our signal are jt(t) and 
jbW ± , bbW^. Also, ti production with both W's decaying leptonically can be a reducible 
background if one of the charged leptons is not detected. This background should be included 
in our analyses when the transverse momentum and pseudo-rapidity of the lepton are in the 
range pr(^) < 10 GeV and |?7(^)| > 2.5. The SM backgrounds always come with W leptonic 
decays with missing neutrinos. To suppress backgrounds, we veto the events with small 
transverse mass of the lepton and missing energy M T = \J {Ert + $t) 2 — {pti + ftr) 2 < 90 



GeV [36|. Furthermore, we take the 6-quark tagging efficiency as 50% and a probability of 
0.4%(10%) for a light (c-quark) jet to be mis-identified as a b jet [34J, |35|. In Fig. @]we show 
the total ti production cross section, with ti — > cb£~(cb£ + ) versus the D mass after 
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FIG. 4: Production cross section of pp — > tiX with ti — > ^ c b+0T vs - mass after basic cuts and 
Mr cut at 14 TeV LHC. Branching fractions for top quark FCNC decay 2x BR t ^ c DD{l— B Rt^cDD) 
are not included while W leptonic decay rate is included. 



cx(pb) 


signals 


jt(i) 


bt(bt) 


jjW ± 


jcW ± + ccW ± 


jbW ± + bbW ± 


ti 


basic cuts 


72 


7.5 


0.32 


2.8 


2.4 


12.7 


0.1 


all cuts 


44 


0.03 


1.6 x 10~ 3 


8.6 x 10~ 3 


0.01 


0.05 


0.05 



TABLE I: ti production cross section with ti — > cb£~ (cbi + ) + $t after basic cuts and Mt cut, 
assuming mo = 20 GeV. Branching fractions for top quark FCNC decay 2 x BRt^cDD^ — 
BRt^cDD) are not included, while the W leptonic decay rate is included. For comparison, the 
background processes are also included with the sequential cut as indicated. 

basic cuts and Mt cut. Assuming mo = 20 GeV, we list in Table U the cross section values 
of our signal and SM backgrounds with basic cuts and Mt cut separately at the 14 TeV 
LHC. One can see that the backgrounds are substantially suppressed. 

After including the appropriate branching fractions for the individual FCNC top quark 
decay, the expected number of events that we are interested in is given by 

N = L x a(pp -> ttX) x 2 x - x BR t ^ cDD x (1 - BR t ^ cDD ), (25) 
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FIG. 5: The 5a discovery limit for BR(t — > cDD) through pp — > ttX with tt -)• £ c b + Xt in the 
BR-m^ plane at 14 TeV LHC with integrated luminosity of 10 fb" 1 (solid) and 100 fb -1 (dashed), 
including all the judicious cuts described in the early section. 

where L is the integrated luminosity. In Fig. [5] we show the 5a signal significance obtained 
in terms of Gaussian statistics, given by the ratio S/ \fB of signal to background events with 
luminosities of 10 fb -1 and 100 fb -1 . Assuming 10 fb -1 luminosity and BR(t — > cDD) > 
10 -4 at 14 TeV LHC, we can expect to observe more than 80 events for mp ^ 60 GeV after 
including all selection cuts and detector effects. With an integrated luminosity of 10 (100) 
fb -1 and the same D mass range, one can explore branching ratios of t — > cDD as low as 
2 x 10" 4 (7 x 10~ 5 ) at 14 TeV LHC. 

V. CONCLUSION 

A stable SM singlet real scalar field, called the D boson, provides a plausible cold dark 
matter candidate that is compatible with the relic abundance measurements. We implement 
this scenario in a two Higgs doublet model (type III) extension which contains tree level flavor 
changing decay t — > cDD mediated by the lightest SM-like Higgs boson h°. The existence of 
D can be explored at the LHC through this FCNC top quark decay, with a branching ratio 
which can approach 10 -3 for m/jO < 2Mw,z- In tt production with tt — )■ cb £~(c b l + ) +J5t, 
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with mjj < 60 GeV and an integrated luminosity of 10 (100) fb 1 at the 14 TeV LHC, one 
can reach 5a significance with a branching ratio BR(t — > cDD) > 2 x 10 -4 (7 x 10 -5 ). 
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